Hydrogen atom in a spherical well: linear approximation 
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We discuss the boundary effects on a quantum system by examining the problem of a hydrogen 
atom in a spherical well. By using an approximation method which is linear in energy we calculate 
the boundary corrections to the ground-state energy and wave function. We obtain the asymptotic 
dependence of the ground-state energy on the radius of the well. 



I. INTRODUCTION 

The problem of a hydrogen atom confined in a sphere 
has quite a long history in quantum physics. It was first 
investigated more than sixty years ago by Michels, de 
Boer, and Bijl Q in their study of the effects of pressure 
on the polarizability of atoms and molecules. This prob- 
lem was then taken up by Sommerfeld and Welker who 
studied the problem in detail and calculated the critical 
radius for which the binding energy becomes zero. Over 
the years there has been a steady flow of papers on this 
and other closely related problems. |3]-|ll|| The model has 
often been used as a test problem for various perturbation 
methods. Using their boundary perturbation method, 
Hull and Julius jl2j obtained a formula which expresses 
the change of energy for the eigenstates in the confined 
system in terms of the corresponding wave functions in 
free space. This method has been improved and gener- 
alized by many authors. |l|,[14| Some variational meth- 
ods have also been used to study this problem. 15 17]] 
Froman, Yngve, and Froman have developed the phase- 
integral method as a general method to attack the prob- 
lem of confined quantum systems and their 1987 paper 
p8[ provides 80 references on this problem. 

In recent years there has been some renewed interest 
on this problem. [ p~9|j20| This is partly driven by the tech- 
nological advances, such as in the field of semiconductor 
quantum dots, (2l]] that have enabled the construction of 
interesting nanostructures which contain a small and con- 
trollable number (1-1000) of electrons. The computation 
of the electronic structure of such systems necessarily has 
to take into account the presence of the finite confining 
boundaries and their influence on the system. 

In this paper we shall study the boundary corrections 
for a hydrogen atom in a spherical well using an approx- 
imation method which is linear in energy. This is a well- 
known method in solid-state physics and has been widely 
used in electronic structure calculations, under the name 
of Linear Muffin-Tin Orbital (LMTO) method, |22l si nce 
its initial introduction by O. K. Andersen in 1975. [ ^3|j24| 
The method is best applied to the calculations of the 
wave functions of a hamiltonian with energies which are 
in close vicinity of the energy of a known wave function. 
To the best of our knowledge, this simple method has 



never been applied previously to this problem of hydro- 
gen atom in an impenetrable sphere. The present paper 
seeks to serve two purposes. First, it presents a new 
approach, which has some pedagogical simplicity, to the 
confined hydrogen atom problem. Second, it offers an an- 
alytically tractable problem from which one can hopefully 
gain some insights into the workings and the accuracy of 
the LMTO method. 



II. LINEAR METHOD 

In this paper we will examine the boundary corrections 
for a hydrogen atom situated at the center of a spherical 
cavity of radius S as shown in Fig.[l| We will assume the 
wall of the cavity to be impenetrable and consider the 
following spherically-symmetric potential: 



V(r) 



-e 2 /r, r < S, 
oo, r > S. 



(1) 



The radius of the cavity will be assumed to be much 
larger than the Bohr radius: S 3> a . In the remainder 
of the paper we shall use the atomic units: 



h = — = 2m = 1. 
2 



The unit of length is the Bohr radius a,Q = Ti 2 /me 2 



the unit of energy is the Rydberg: Ry 



V2ao 



(2) 

and 
13.6 



eV. The Schrodinger equation takes the following form: 



£T*(r) 



*(r) = E$(r) 



(3) 



The wave function ^(r) satisfies the Schrodinger equa- 
tion for the hydrogen atom for r < S, in particular it 
should still be regular at the origin. The only difference 
from the free-space case is that now we have to impose a 
different boundary condition: the wave function should 
vanish at r = S instead of at r = oo . 

For S ^> do, the changes in the ground-state wave 
function and energy due to the presence of the wall are 
expected to be "small" because the wave function is con- 
centrated at the center of the cavity, far away from the 
confining wall. 
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FIG. 1. Hydrogen atom in a spherical well of radius S. 

In free space, i.e. in the absence of the confining cavity, 
the hydrogen atom has the familiar Rydberg spectrum: 



-~2, n = l,2, 



In the presence of the cavity, we write 
E n — e n + Ae„ . 



(4) 



(5) 



We use small letters (e, ip, etc.) to denote quantities for 
the free-space problem and capital letters (E, "J, etc.) 
for the corresponding quantities in the cavity problem. 
The dimensionless parameter (Ae„/e„) is expected to be 
small for n 2 ao <C S. In the linear method, the (unnor- 
malized) wave function at energy E n is approximated by 



*(E„,r) = ip(e n , r) + Ae„ ip(e n , r). 



(6) 



Here ip(£ n ,r) is the derivative with respect to energy of 
ip(e, r) evaluated at e = e n : 



[d^(e,v)/de] (e 



(7) 



The eigenfunctions in the cavity problem are then ob- 
tained by imposing the boundary condition at r = S: 

S,r) = 0, (8) 

which gives an expression for the energy correction: 



ip(s„,S,r) 



(9) 



Here r = (9, </>) is a unit vector in the direction of r. 

To apply the linear approximation method we need the 
general solution to the Schrodingcr equation at an arbi- 
trary energy E. Since we are dealing with a spherically- 
symmetric system, we can separate the variables: 



V(r) = R(r)Y lm (r). 



(10) 



The resulting radial differential equation is 

1(1 + 1) 



d 2 R 
dr 2 



2dR 

r dr 



2 

E+- 

r 



R = 0. 



Transforming the variables by defining 
lu = y—E, p = 2u>r, 



(11) 



(12) 



and using the following trial functional form |2q | 

R(p) = P l e-el 2 u(p), (13) 



then gives us the following differential equation 25 

1 



pu 



2(1 + 1)- p 



1 + 1- 



u = 0, (14) 



which is the equation for the confluent hypergeometric 
function. The general solution of this equation, which is 
regular at the origin, is p5fl 



u(p)=A 1 F 1 (l + l-^-;2l + 2- p y 



(15) 



where A is a normalization constant. The radial part of 
the general solution to the Schrodinger equation Eq.(|^) 
with energy E = —lu 2 therefore is 

Rl(w, r)= A (2ior) l e^ r 1 F 1 (l + 1 ; 21 + 2; 2cor 



(16) 

The free-space solution is obtained by requiring that 
R(r) — > as r — > oo. From the properties of the hyperge- 
ometric functions, j25) this can only happen if (1+1 — 1/uS) 
is a negative integer or zero. This implies that 



1 



-=n, I = 0,1, ...,(n- 1), 

CO 



(17) 



with n a positive integer. This directly leads to the Ryd- 
berg spectrum in Eq.(^). 

The function Ri(w,r) is plotted in Fig.|| for I — and 
u> = 1, 0.98, and 0.50. The u> = 1 curve is the ground- 
state wave function of the hydrogen atom in free space 
and is nodeless. Here a node of Ri(u,r) is defined to 
be a value of the argument r which gives zero value for 
the function Ri(u>, r). As w is reduced below 1, the wave 
function acquires a single node which moves from r = oo 
to r — 2ao at lu — 0.50. where it becomes the (n, I) — 
(2, 0) eigenstate of the hydrogen atom in free space. One 
therefore can obtain the ground-state wave function and 
energy of the hydrogen atom in a cavity of radius S by 
numerically searching for the energy which gives a wave 
function with a single node at r = S. This provides a 
useful comparison for our approximation. 

Since the spherical harmonics are independent of the 
energy we can recast Eq.(||) into 



2 
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FIG. 2. The function rRi(uo,r) as a function of r/ao for 
I = and u = 1, 0.98, and 0.50. The w = 1 curve is nodeless. 
As w is decreased from 1 to 0.50, the node of the wave function 
moves from r — oo to r = 2ao. 



A£ n ; — 2w n — : — - — -. 

Ri(u n ,S) 



where u> n = \J—e n and 

Rl(cJ n ,S) = [8Ri(u;,S)/du J } (w = w„). 



(18) 



(19) 



Substituting the radial function Ri(u,r) in Eq. (|16[) into 
Eq.(|l8|) then gives us an explicit formal expression for 
Ae„ which should be valid for R 3> n 2 a$. Note that 
the presence of the finite boundary lifts the azimuthal 
degeneracy of the states with different orbital quantum 
number I (and the same radial quantum number n). J2(| 
As in the case of the screened Coulomb potential, this 
occurs because one no longer deal with the pure Coulomb 
potential. [f^H 

To gain an insight into Eqs.(|lq)-(|l9|), we shall consider 
the ground state (n — 1), which is a special case of the 
zero angular momentum (Z = 0) states. We have 



1 



iFi(l-~;2;2wr 



Rq(u, r) = A e 
For the ground state (n — 1), this is 

i?o(l, r) = A e- r iFi (o; 2; 2rj = A e 



(20) 



(21) 



We are interested in obtaining a simple analytical ex- 
pression of the correction to the ground-state energy for 
S 3> do, therefore we need to calculate the limiting form 
of R{){uj, r) for r ^> oq. The asymptotic expansion of the 
hypergeometric function i-Fi(a, b, z) for large z is [E9[ 



iFi{a,b,z) 
T(b) 



h(a,b,z) | cZ „q-6 h(a,b,z) 



r(6- 



with 



h(a,b,z) 



y-J (a)n(l + « - fr)r 

n=0 



n 



r(a) 



+ 0(\z\- R ), (23) 
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FIG. 3. Dependence of the ground-state energy of a hy- 
drogen atom confined in a spherical cavity on the radius of 
the cavity S. The topmost curve is the exact result which is 
obtained by numerically searching for the node of the wave 
function. The middle curve is obtained from the linear approx- 
imation, Eq.(|l^), using the exact wave function Eq. (pl|) . The 
lowest curve is obtained using the limiting formula Eq. ( |3l| ) . 



I 2 (a, b, z) = I! {b - ^ - a)n -L + 0[\z\~ R ). (24) 

n=0 

The Pochhammer symbol (a)„ is defined by (25) 

T(a + n) 



(a)„ = a(a + 1) • • • (a + n — 1) = 



(25) 



We need to calculate the derivative of this function at 
a = — with u> = 1. In this case the dominant term 
comes from the derivative of T(a) in the second term in 
Eq. (p2[) . The first term can be neglected because it does 
not have the exponential term e z which dominates the 
derivative at large distances. Keeping only the largest 
term, we get 



d_ 

Oa 



iFi(a,b,z) 



-e z z a - b 



T(b)I 2 (a,b,z) (26) 



Here ip{ a ) is th e digamma function: ip(a) = T'(a)/T(a) 
p9[ Its ratio with T(a) as a — > is 



, ^(«) , 
lim — 7-t = lini 



-7 - 1/a 



= -1, 



a^o r(a) a^o — 7+l/a 
where 7 is the Euler constant. This then gives 
d 



(27) 



da 



Fi(a,b,z) (a^0)^e z z a - b T{b)I 2 (a,b,z). (28) 



3 



Using this expression, and keeping only the first two 
terms in l2(a,b,z), we can obtain the limiting form of 
Rq(u>, r) at large r and u> — > 1: 



£o(w,r) 



Ae 



oj 2 I (2wr) 1 + 1 ^ 



1 



r(2 + i/uj) 

2wr r(l/w) 



Exactly at cj = 1 , this expression becomes 



ifc(l,r) 



4r 2 



(29) 



(30) 



Finally, using this equation and Eq.(|2l|) in Eq. (|Ls|) , we 
get the boundary correction to the ground-state energy: 



Ae (S) S3 8S(S-l)e 



-2S 



S > CLq. 



(31) 



III. DISCUSSION 



Fig.^ displays the asymptotic dependence of the energy 
correction on the radius of cavity, Eq. ( |3l| ) , together with 
the exact curve and the curve obtained from the linear 
approximation method, Eq.(|l8|), using the exact wave 
function Eq.(21). It is seen that the asymptotic formula, 
Eq. (|3l|) , is fairly accurate for radii greater than about 
four Bohr radius. Note that the exact energy at S = 2ao 
is equal to \ Ry, which is the energy of the first excited 
state (n,l) = (2,0) of the hydrogen atom in free space. 
This is because the corresponding wave function has a 
node at r = 2a as can be seen in Fig.^. 

The asymptotic formula Eq.(|3l"|), which is the limit 
curve in Fig.|J is a "double-approximation" to the ex- 
act curve. It is an asymptotic form of the linear curve, 
Eq. (|l8|) , valid for large values of S/oq. The linear curve 
itself is an approximation, linear in energy, to the ex- 
act curve. For small values of S/clq, and within the lin- 
ear approximation method, one has to use Eq.(|lg|) which 
in general, unfortunately, does not correspond to a sim- 
ple analytic expression. This does not pose a problem 
in actual electronic-structure applications because there 
the wave function and its energy derivative are computed 
numerically. In this paper, for pedagogical purposes, we 
have calculated the asymptotic formula, Eq.(|3l|), which 
does correspond to a simple analytic expression. 

Knowing the dependence of the ground-state energy 
on the cavity radius, Eq.(|3l]), allows us to calculate the 
pressure needed to "compress" a hydrogen atom in its 
ground state to a certain size. This is given by 



P(S) 



dAs _ 4e 
dV ■ 



-2S 



(>-§) 



(32) 



At S = 4a this has a value of 2.13 x 10" 4 eV/ajj = 
1.47 x 10 4 GPa. At this radius, the change of the ground- 
state energy is 0.032 Ry which is only three percent of 
the binding energy of a free hydrogen atom. 



The information on the effects of the boundary on the 
wave function of the atom can also be used to study 
the influence of the boundary on other properties of the 
atom, e.g., the spin-orbit coupling energy. It is also inter- 
esting to calculate the changes in the wave function and 
energy of the atom when it is displaced from the center of 
the cavity, and the force that will push it back towards 
the center. The linear method also seems to be well- 
suited for the analysis of the "soft-cavity" case where we 
have a finite potential outside the cavity, instead of the 
infinite potential considered in this paper. These topics 
will be examined in future works. 

In conclusion, we have used a linear approximation 
method to calculate the asymptotic dependence of the 
ground-state energy of a hydrogen atom confined to a 
spherical cavity on the radius of the cavity. The bound- 
ary correction to the energies of the excited states can be 
obtained using the same method. 
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